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pH Here I review measurements of V2, the second component in a Fourier decomposition 

of the azimuthal dependence of particle production relative to the reaction plane in 
heavy-ion collisions. V2 is an observable central to the interpretation of the subsequent 
expansion of heavy-ion collisions. Its large value indicates significant space-momentum 
^ ' correlations, consistent with the rapid expansion of a strongly interacting Quark Gluon 

nJ , Plasma. Data is reviewed for collision energies from y^sTTjv = 2 to 200 GeV. Scaling 

1 observations and comparisons to hydrodynamic models are discussed. 

^ ■ 

p ; 

ly^ i 1. Introduction 

o 



Collisions of heavy nuclei have been exploited for decades to s earch for and study 
. the transition of hadronic matter to quark gluon plasmaPJ^I Jn these collisions, 

the extended overlap area, where the nuclei intersect and initial interactions occur, 
does not possess sphrerical symmetry in the transverse plane. Rather, for non- 
central collisions, the overlap area is roughly elliptic in shape. If individual nucleon- 



X 

' nucleon collisions within the interaction region are independent of each other {e.g. 

point-like) and no subsequent interactions occur, this spatial anisotropy will not be 
reflected in the momentum distribution of particles emitted from the interaction 
region. On the other hand, if the initial interactions are not independent, or if there 
are subsequent interactions after the initial collisions, then the spatial anisotropy 
can be converted into an anisotropy in momentum-space. The extent to which 
this conversion takes place allows one to study how the system created in the 
collision of heavy nuclei deviates from a point-like, non-interacting system. The 
existence and nature of space-momentum correlations is therefore an interesting 
subject in the study of heavy ion collisions and the nature of the matter created 
in those collisions''^. Fig. [T] shows an illustration of the possible stages of a heavy- 
ion collisions starting with some initial energy density deposited at mid-rapidity, 
followed by a QGP expansion, a hadronization phase boundary, a kinetic freeze-out 
boundary and finally the observation of particle trajectories in a detector. 
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Nuclear collisions and tlie QGP expansion 




Fig. 1. A schematic diagram of the expansion after an ultra-relativistic heavy-ion coUision. 



One can consider a number of ways to study space-momentum correlations: e.g. 
two-particle correlations-' and HBTSi^. In this review we discuss elliptic flow V2 ; an 
observable that has been central in the interpretation of heavy-ion data and QGP 
formation . Given the predominantly elliptic shape of the initial overlap region, it 
is natural to ask whether this shape also shows up in the distribution of particles 
in momentum-space. Fig. [5] shows a schematic illustration of the conversion of 
coordinate-space anisotropy to anisotropy in momentum-space. The left panel shows 
the position of nucleons in two colliding nuclei at the moment of impact. The overlap 
region is outlined and shaded. A Fourier decomposition can be used to describe the 
azimuthal dependence of the final triple momentum-space distribution^ 

(fN 1 dN 

(1) 

where 4> is the azimuth angle of the particle, y the longitudinal rapidity variable, 
Pt the transverse momentum, and ^' is the reaction plane angle defined by the 
vector connecting the centers of the two colliding nuclei. Positive V2 implies that 
more particles are emitted along the short axis of the overlap region. To study 
the extent to which space-momentum correlations develop in heavy-ion collisions, 
one can measure the second component V2 and compare it to the initial spatial 
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Fig. 2. Schematic illustrations of a y^s^Tjv ~ '^^^ GeV Au+Au collision with a 6 fm impact 
parameter. The left panel shows the nucleons of the two colliding nuclei with an ellipse outlining 
the approximate interaction region. The right panel shows a momentum-space representation of V2- 
The average radius of each successive ring represents the px of the particles while the anisotropy of 
the ring represents the magnitude of V2 ■ The highest pj- particles (outer-ring) exhibit the strongest 
V2 while the lowest px particles (inner-ring) exhibit a vanishingly small V2. 

eccentricitjEHni The right panel of Fig. [5] shows the final azimuthal distribution 
of particles in momentum-space. The curves represent the anisotropy at different 



2^2 (pt) c,os{atan2{py , Px))) ■ The goal of V2 measurements is to study how the initial 
spatial anisotropy in the left panel is converted to the momentum-space anisotropy 
in the right panel. In this review, a summary of V2 data for different colliding 
systems, different center-of-mass energies, and different centralities is given. 

This review will focus on results from the first four years of operation of the 
Relativistic Heavy Ion Collider (RHIC). We start with a brief discussion of the beam 
energy dependence of V2 and some ideas about what physics might be relevant. Even 
before considering physics scenarios to explain how a space-momentum correlation 
develops, one can see that to interpret V2 it is important to understand the initial 
geometry and how it varies with the collision centrality and system-size. Since, the 
concept of the reaction-plane is so central to the definition of V2 and eccentricity is so 
central to it's interpretation, I discuss the two in a sub-section below. Then a review 
of RHIC data is provided. This will include the dependence of V2 on center-of-mass 
energy, centrality, colliding system, pseudo-rapidity, p^, particle mass, constituent 
quark number and various scaling laws. In the following section, I will discuss 
comparisons to models and the emergence of the hydrodynamic paradigm at RHIC. 
Particular emphasis is given to uncertainties in the model comparisons. In that 
section I will also discuss current attempts to extract viscosity and future directions 
of investigation. 

Voloshin, Poskanzer, and Snellings recently wrote a review articl^^on collective 



Pt values measured 
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phenomena in non-central nuclear collisions that deals with a similar subject matter. 
That article provides valuable detail on technical aspects of measuring V2- In this 
review I will attempt to avoid duplicating that work by discussing interpretations 
of V2 more extensively and refer the reader to that review where appropriate. 



1.1. Two Decades in Time and Five Decades in Beam Energy 
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Fig. 3. The beam energy dependence of elliptic flow measurements. The RHIC and E877 data are 
for charged hadrons, independent of species. The NA49 data is for charged pions. The E895 data 
are for protons and the FOPI data are for atomic number Z=l. At each energy, the sample of 
particles is close to the total charge. Positive V2 values indicate that particles tend to be more 
aligned with the reaction-plane (in-plane). RHIC and SPS data suggest a smooth trend of in-plane 
V2 growing with log(^s„„) above Eteam/A ~ 20 GeV or y^sTTiv" ~ 6 GeV. 



Positive values of V2 imply that particles tend to be produced more abundantly 
in the x— direction than in the y— direction. This is referred to as in-plane flow. 
Fig. [3] sh ows t'2 measured in an interv al of beam energies covering five orders of 
magnitud J^'^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l For Ebeam/A ranging from approximately 
0.12 — 5 GeV (1.92 < y/s^ < 3.3 GeV), V2 is negative. For this energy range, spec- 
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tator protons and neutrons are still passing the interaction region while particles are 
being produced. Their presence inhibits particle emission in the in-plane direction 
leading to the phenomenon termed squeeze-out. At still lower energies, V2 is posi- 
tive as the rotation of the matter leads to fragments being emmitted in-plane. At 
this energy beam rapidity and mid-rapidity are essentially indistinguishable with 

Vbeam < 0.41 UnitS. 

In-plane flow: As the beam energy is increased, the nuclei become more 
Lorentz contracted and the time it takes the spectators to pass each other decreases. 
It was predicted by OUitraultP that at high enough beam energy, the squeeze-out 
phenomena would cease and V2 would take on positive values. Positive V2 values 
were measured at the AGS for energies above Ebeam/A = 5 GeV {y/s^^ = 3.3 GeV). 
For energies above Eheam/^ ~ 20 GeV (y^s^ = 6.3 GeV), V2 exhibits a steady log- 
linear increase: V2 w 0.01 + 0.0042 log(i;fceamM) or V2 « 0.008 -I- 0.0084 log(^^) 
where the data represented are from intermediate impact parameter A-l-A collsions. 
It appears therefore that RHIC V2 data may be part of a smooth trend that began 
at SPS energies. This trend was noted previously at least oncJ23l Understanding 
the physics that underlies that trend is one of the challenges of heavy-ion physics. 

One class of models that has provided an illustrative reference for heavy-ion col- 
lisions are hydrodynamic models which are used to .^^gdel the expans m atter 
remaining in the fireball after the initial collision!p^ ^^ l ^^ l ^'' l ^^ l ^^ l ^" l ^l | ^^ l '^^ l. This 
model can be used to determine how matter with a vanishingly small mean free 
path would convert the initial eccentricity into V2- These models typically treat all 
elliptic flow as arising f rom the final state expansion rather than from some initial 
state 

,ffe,y34|35|36|37| 

In the hydrodynamic models, large pressure gradients in 
the in-plane direction lead to a preferential flow of matter in the in-plane direc- 
tion. In this review, we will use hydrodynamic models as a convenient reference. 
Other models providing a valuable refer ence for measurements inc lude hadronic and 
partonic cascades and transport modelj^S|41| 42 | 43 | 44 | 39 | 40 | 45 | 46 | Additionally, the 
blast-wave model provides a successful parametrization of low pr heavy -ion d ata, 
including V2, HBT, and spectra in terms of several freeze-out parameterJ^OIlZl, 



1.2. Initial Geometry: The Reaction Plane and Eccentricity 

In the collision of two symmetric nuclei, a unique vector (the y-axis) can be defined 
by applying the right-hand-rule to the momentum vector of one nucleus and the 
vector pointing to the center of the other nucleus. The y-axis is a pseudovector. The 
reaction-plane is then the plane perpendicular to the y-axis containing the points 
at the center of the two nuclei. The reaction-plane and the right-handed coordinate 
system are illustrated in Fig. ID The figure contains perspective illustrations of two 
nuclei approaching with an impact parameter of 6 fm. The impact parameter is 
the distance between the centers of the two nuclei at the moment of their closest 
approach. The two nuclei in this illustration are Lorentz contracted by a Lorentz 
gamma factor of 10 which roughly corresponds to the appropriate gamma for top 
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Fig. 4. Schematic illustrations of a y/s~^^ = 200 GeV Au+Au collision with a 6 fm impact 
parameter. The right-handed coordinate systems defined by the momentum of the nucleus (here 
the z-axis) crossed with the vector pointing to the center of the approaching nucleus (the x-axis) 
is shown for each nucleus. The reaction-plane is the plane normal to the y-axis, containing the 
centers of the two colliding nuclei. The nucleons are distributed inside the nucleus according to a 
Woods-Saxon distribution. The nuclei are Lorentz contracted in the z-direction. 



SPS energies. 

The reaction-plane is not directly observed in experiments, however, and this 
introduces a systematic uncertainty into the measurement of V2. One often relies 
instead on indirect observations to estimate t;248 49 50 51 p^j, example, when 
forming two particle azimuthal correlations such as ^(^3^1 ^ non-zero V2 value 
will lead to a modulation in Acf) ^ (f>i — <j>2 of the form 1 + 2{v2) cos(2A(/)). Fig. \5\ 
shows the correlation function for hadrons produced at mid-rapidity at RHIO^^IMI^ 
The panels show different centralities. The area normalized correlation function is 

where Ygl^^{A(j)) and Y^^^^{A(I)) are, respectively, the uncorrected yields of pairs 
in the same and in mixed events within each data sample. C(A(/)) shows a clear 
cos(2A(/)) dependence. We note here that what is measured in these correlation 
functions is — + <^l., in anticipation of a discussion of V2 fluctuations. 

V2 will not be the only contribution to the azimuthal dependence of the two- 
particle azimuthal correlations. Other processes that are not related to the reaction- 
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Fig. 5. Charged hadron correlation functions in Au+Au collisions at ^fs^ 
centrality intervals and two ranges. 



200 GeV for two 



plane can give rise to structures in the shape of the two-particle At/) distribution as 
well. These non- reaction plane contributions are commonly called " non-flow" . The 
subject of non-flow is an important one and will be discussed throughout this review. 
The contribution of non-flow can be seen more clearly by looking at very peripheral 
collisions or by selecting high momentum particles to increase the chance that a 
particular pair of hadrons are correlated to a hard scattered parton (jet). Fig. [6] 
shows the correlation function for higher momentum particlei^. The solid line 
shows what a pure V2 correlation would look likc^^. The difference bet ween thos e 
curves and the data are often taken as a measurement of jet correlationiSSl57l58]_ 

Even if the reaction-plane were known with precision, there is no first principles 
calculation of the initial matter distribution in the overlap region, so the eccentricity 
is uncertain. Various models can be used to calculate the initial spatial eccentricity 
which can then be compared to V2. Defining the y— axis according to the right- 
hand-rule, the eccentricity Es is traditionally calculated as: 

where the average represents a weighted mean. Other eccentricity definitions have 
also been considered^. The weights can be some physical quantity in a model such 
as energy or entropy density, or simply the position of nucleons participating in 
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Fig. 6. Charged hadron correlations for a variety of centrality intervals. The correlation function 
is formed between two samples of hadrons based on their p^: the "trigger" particle sample is 
selected with 2.5 GeV/c < pT < 4.0 GeV/c while the "associated" particle sample is selected with 
1 GeV/c <pt < 2.5 GeV/c. 



the collision. One popular method for calculating the eccentricity is to use a Monte 
Carlo Glauber model. Details can be found in a recent RevieWSQj, jn that model, a 
finite number of nucleons are distributed in a nucleus according to a Woods-Saxon 
distribution. Then two nuclei are overlaid with a fixed impact parameter and the 
X and y positions of the participating nucleons is determined based on whether 
the nucleons overlap in the transverse plane; each nucleon is considered to be a 
disk with an area determined by the ^/s dependent nucleon-nucleon cross-section. 
The X and y coordinates of the participating nucleons are then used to calculate 
the eccentricity. Those nucleons that do not participate in this initial interaction 
are called spectators. One can anticipate that due to the finite number of nucleons 
in this model, the initial geometry will fluctuate. Ot her mod els used to determine 
the initial matter distributioii including HIJINgISII, NEXU^, and Color Glass 
Condensate also reach the same conclusions; the initial overlap 

region is expected to be lumpy rather than smooth. Fig. [7] shows the gluon density 
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in the transverse plane which is probed by a .2 fm quark-antiquark dipole at two 
different x values in the IPsat CGC mode^^ {x — 2pt/v^s^^ is 10^^ in the left 
panel and 10^"^ in the right panel). The lumpiness is immediately apparent. 



Fig. 7. Gluon density in the transverse plane when the nucleus is probed at different x values by 
a 0.2 fm quark-antiquark dipole in the IPsat CGC model. 



Until recently howeveJ ^'^ l ^^ i V2 data was compared almost exclusively to cal- 
culations assuming an infinitely smooth initial matter distribution (for example in 
initializing a hydrodynamic expansion). Improving on that approximation ma y be 
important for understanding the shape expected for the ^^^^^ distribution^ 
This distribution is also investigated in heavy ion collisions in order to search for 
jets. In any scenario where space-momentum correlations develop, the correlations 
and fluctuations in the initial geometry can be manifested in the distribution 
and understanding these correlations is important for interpreting heavy-ion colli- 
sions. Fluctuations in the initial geometry have also led to the idea of measuring 
parti cle dis tributions relative to the participant-plane rather than the reaction- 
planPEi. The participant-plane is defined by the major axis of the eccentricity 
which, due to fluctuations, can deviate from the reaction-plane. The eccentricity 
relative to the participant-plane is a positive definite quantity and is always larger 
than the eccentricity relative to the reaction-plane; the participant-plane is defined 
by rotating to the axis that maximizes the eccentricity. Fig. [8] shows the event- 
by-event distribution of the standard eccentricity (left panel) and the participant 
eccentricity (right panel) as a function of impact parameter determined from a 
Monte-Carlo Glauber calculation. The fluctuations in this model are large as illus- 
trated by the widths of the distributions. The relationship between the different 
definitions of eccentricity and their fluctuations are explained clearly in two recent 
paperPEB. 
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Fig. 8. The distribution of eccentricity on an event-by-event basis when calculated relative to the 
reaction-plane (left) and the participant plane (right). 

Different models for the initial matter distribution yield different estimates of (e) 
and ^/J^^. The deviations in the (e) for different models can be of the order of 30% 
and strongly centrality dependent!^. That uncertainty in (e) leads to an inherent 
uncertainty when comparing models to V2/£- This level of uncertainty becomes 
important when attempting to estimate transport properties of the matter based 
on comparisons of the observed V2 to the initial e. 



2. Review of Recent Data 

The first paper published on RHIC data was on elliptic flow in ^/Tmn = 130 GeV 
Au+Au collisions Fig. [9] shows that data on the centrality dependence of V2. 
The values of V2 reach a maximum of approximately 6% for peripheral collisions 
where the initial eccentricity of the system is largest. That value is 50% larger 
than the values reached at SPS energieJ^H and the V2 values are a factor of two 
larger than the those predicted by the RQMD transport-cascade modeS^. For cen- 
tral collisions, the measurements approach the zero mean-free-path limit estimated 
from the eccentricity shown in the figure as open boxes. The boxes represent the 
eccentricity scaled by 0.19 (bottom edge of the boxes) and 0.25 (top edge of the 
boxes). Those values are chosen to represent the typical conversion of eccentricity 
to V2 in hydrodynamic models. At lower ^ysNN energies, the RQMD model pro- 
vided a better description of the data, while hydrodynamic models significantly 
over-predicted the data. The conclusion based on this early comparison, therefore, 
was that heavy-ion collisions approximately satisfy the assumptions made in the 
hydrodynamic models: 1 ) ze ro mean-free-path between interactions, and 2) early 
local thermal equilibriunl231. These conclusions remain at the center of scientific 
debate in the heavy-ion community. 

In Fig. [ini V2 is scaled by model calculations of the initial eccentricity and 
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Fig. 9. First measurements of tJ2 versus centrality at RHIC. Positive values are observed which are 
largest for events with the largest eccentricity and decrease for more central, symmetric collisions. 
The trend with centrality clearly indicates a space-momentum correlation driven by the eccentric- 
ity of the initial overlap zone. Centrality is expressed in terms of the observed multiplicity of a give 
event relative to the highest multiplicity observed n^i^/nmax- Data are compared to eccentricity 
scaled by 0.19 (bottom edge of the boxes) and 0.25 (top edge of the boxes). The values are chosen 
to represent the typical conversion of eccentricity to V2 in a hydrodynamic model. 



plotted versus transverse particle density ^ ■ This facilitates comparisons of 
V2 across different ^s„„ energies, collision centralities and system-sizeJ^Hini, 

For the case of ballistic expansion of the system - that is an expansion for 
which the produced particles escape the initial overlap zone without interactions - 
V2 should only reflect the space-momentum correlations that arise from the initial 
co ndit ions. Those can exist in the case that the initial interactions are not point- 
likd^ but rather involve cross-talk between different N + N interactions within the 
overlap zone. The opposite extreme from the ballistic expansion limit is the zero 
mean-free-path limit represented by ideal hydrodynamic models. Lacking a length 
scale, the zero mean-free-path models should not depend on system-size and instead 
should be a function of density. 

The measurements of V2 are expected to rise from values near the ballistic 
expansion limit and asymptotically approach the zero mean-free-path limit as the 
density of the system is increased. Data in Fig. \W\ exhibit such a behavior with 
the most central collisions at full RHIC energy apparently becoming consistent 
with the hydrodynamic model. This conclusion however depends on the model 
calculations for the initial eccentricity and on the assumption that the observed V2 
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Fig. 10. V2 scaled by a Monte Carlo Glauber model calculation of the initial overlap eccentricity. 
The ratio is plotted versus transverse particle density {1 / S)dN / dy , where S = tt-^/ (x^) {y'^) is 
a weighted average area calculated with the same model as the eccentricity. Data are taken from 



different 



values and different centralities. Plotted in this format, the data suggest V2/s 



for different energies and overlap geometries is determined by the transverse particle density, and 
approaches a zero mean-free-path hydrodynamic limit for most central top energy collisions at 
RHIC. This conclusion is not universally accepted and is still being investigated. 



dominantly arises from an expansion phase where anisotropic pressure gradients are 
the origin of the space-momentum correlations. Different models for the eccentricity 
yield e results that deviate both in their centrality dependence and in their overall 
magnitude. Reasonable models for the eccentricity can easily give magnitudes 30% 
larger than those used in Fig. [10] with a stronger centrality dependence. The ratio 
V2/S can therefore be smaller than what is shown and have a different shap4^. 
Given this level of uncertainty, the conclusion that heavy-ion collisions at ^/s^ — 
200 GeV approximately satisfy the assumptions made in the hydrodynamic models 
i.e. early local thermal equilibrium and interactions near the zero mean- free path 
limit, would be more convincing if an asymptotic approach to a limiting value were 
observed. Rather, for the eccentricity calculation used in Fig. [TOl the data suggest 
a nearly linear rise with no indication of asymptotic behavior. 

In the hydrodynamic picture, one might also expect that V2/e versus will 
be sensitive to the equation-of-state of the matter formed during the expansion 
phase. Since V2 is expected to reflect space-momentum correlation developed due 
to pressure gradients and ^ is a measure of the transverse particle density, U2 /e 
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Fig. 11. The QCD equation-of-state (pressure over energy density versus the fourth power of the 
energy density) as determined in Lattice calculations. 



versus ^ could be considered as a proxy for the pressure versus energy density 
or the equation-of-state. It's difficult to identify in the data the features that are 
expected in the equation-of-state. Fig. [TT] shows the equation-of-state calculated in 

•JO 

a recent lattice QCD calculation . The onset of the QGP phase is seen to lead to 
an increase in the pressure as the energy density is increased above a critical value. 
The energy density in heavy-ion collisions is often estimated from the Bjorken 
formulaM-^ 

1 dEr ... 
'^= = A-r^ 

which depends only on the experimentally accessible quantity , the overlap area 
of the nuclei and r, the unknown formation time which is often assumed to be 1 
fm. The Bjorken estimate for the energy density is closely related to the transverse 
particle density {I / S){dN / dy) . 



2.1. Differential Elliptic Flow 

In addition to studying how V2 integrated over all particles depends on the central- 
ity or i/s^Tiv of the collision, one can study how V2 depends on the kinematics of 
the produced particles (differential elliptic flow) . Fig. [12] shows the centrality and 
pseudo-rapidity dependence of V2 for 200 GeV Au-f Au collisions^. V2 is largest 
at mid-rapidity where the transverse particle density is largest and then falls off 
at larger jryj values. This behavior is therefore consistent with the trends seen in 
integrated V2 where V2/£ appears to increase with increasing transverse particle 
density. 
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Fig. 12. 1^2 (^) for 200 GeV Au+Au collisions for three centrality intervals. The inset shows the 
ratio of V2{ri) for central collisions over peripheral collisions. Open symbols are data reflected to 
negative rj. 



The fall off oiv2{'ri) with increasing \rj\is common to the three centrality intervals 
studied. The inset of the figure shows the ratio of V2 in peripheral over central 
collisions. Within errors the ratio is flat indicating a similar shape for all centralities 
with V2{ri) only changing by a scale factor. Scaling of W2('7) for different energies 
and system sizes will be discussed in a later section. 

Fig. ll3l sh ows ti2 for a variety of particle species as a function of their transverse 
momentum p ^ | 76 | 77 | 78|79 | 80 |81J the region below pT ^ 2 GeV/c, V2 follows mass 
ordering with heavier particles having smaller V2 at a given px- Above this range, 
the mass ordering is broken and the heavier baryons take on larger V2 values. 

A hydrodynamic model for V2{pt) is also shown which describes the V2 in the 
lower pt region well. This mass ordering is a feature expected for particle emission 
from a boosted source. In the case that particles move with a collective velocity, 
more massive particles will receive a larger px kick. As the particles are shifted to 
higher pT, the lower momentum regions become depopulated with a larger reduction 
in the direction with the largest boost (in- plane). This reduction reduces V2 at a 
given Pt, with the reduction largest for more massive particles. Note that this does 
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Fig. 13. V2{pt) for a variety of identified particle species and inclusive charged hadrons. Data are 
for Au+Au collisions at 200 GeV averaged over all centralities. At low momentum 112 (pt) exhibits 
mass ordering while at larger pT the identified particle V2 appears to be grouped according to 
constituent quark number. The mass ordering at low pT is approximately reproduced by the 
hydrodynamic calculation. 



not imply that the more massive particles have a smaller integrated V2 value, and in 
fa ct the opposite is true. Fig. [Ml shows V2 for identified particles integrated over all 
p^ 82|83 | 'pjjg integration shows that V2 increases with particle mass. This is because 
the more massive particles have a larger (px) and V2 is generally increasing with pT 
in the px region where the bulk of the particles are produced. The hydrodynamic 
model also exhibits this trend. 



2.1.1. Identified Particle V2{pt)'- RHIC versus SPS 

Fig. [13 shows pion and proton V2 from y^s^ ~ 62.4 Au+ AiJSIl and 17.3 GeV 
Pb+Pb colhsionJ^. The centrality intervals have been chosen similarly for the 
17.3 GeV and 62.4 GeV data. The STAR data at 62.4 GeV are measured within 
the pseudo-rapidity interval I77I < 1.0 and the 17.3 GeV data are from the rapidity 
interval < y < 0.7. These intervals represent similar y/ybeam intervals. It has been 
shown that V2 data for pions and kaons at 62.4 GeV are similar to 200 GeV data; 
the 62.4 GeV data only tending to be about 5% smaller than the 200 GeV data. 

Appreciable differences are seen between the 17.3 GeV and 62.4 GeV data. At 
PT > 0.5 GeV/c, for both pions and protons, the V2 values measured at 62.4 GeV 
are approximately 10%-25% larger than those measured at 17.3 GeV. Although the 
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Fig. 14. pt integrated V2 versus mass for a variety of identified particle species. Yellow bands 
show hydrodynamic calculations with different freeze-out temperatures. 



magnitude of V2 is different at the lower energy, the systematics of the particle-type 
dependencies are similar. In particular, pion V2 and proton V2 cross over each other 
at Pt near 1.7 GeV/c for ^s„„ = 17.3, 62.4 and 200 GeV data. Due to the Hmited 
kinematic range covered by the 17.3 GeV data, it's not possible to determine if the 
V2 of baryons at > 2 GeV becomes larger than that for the lighter mesons. 

The increase in the magnitude of V2 from 17.3 GeV to 62.4 GeV and the sim- 
ilarity of 62.4 GeV V2 to 200 GeV V2 has been taken as a possible indication for 
the onset of a limiting behavioJ^. In a coUisional picture, a saturation of V2 could 
indicate that for ^s„„ at and above 62.4 GeV the number of collisions the system 
constituents experience in a given time scale can be considered large and that hy- 
drodynamic equations can therefore be applied. Hydrodynamic model calculations 
of V2 depend on the model initialization and the poorly understood freeze-out as- 
sumptions. As such, rather than comparing the predicted and measured values at 
one energy, the most convincing way to demonstrate that a hydrodynamic limit 
has been reached may be to observe the onset of limiting behavior with ^s„„ . For 
this reason, V2 measurements at a variety of center-of-mass energies are of interest. 
Fig. [T5l shows that when the 17.3 and 62.4 GeV V2{pt) data are compared within 
similar |j/|/j/6eam intervals, the differences between V2{pt) within the data sets may 
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Fig. 15. V2{pt) for pions and protons at 



62.4 and 17.3 GeV. 



be as small as 10%-15%. As such, a large fraction of the deviation between the 
SPS data and hydrodynamic models arises due to the wide rapidity range covered 
by those measurements {v2 approaches zero as beam rapidity is approached^^ , 
increased (pr) values at RHIC and the larger V2 values predicted for the lower 
colliding energy by hydrodynamic models. 

2.2. High px 

At higher px, ^2 no longer rises with px and the mass ordering is broken. Above 
Pt ^ 2 GeV/c the more massive baryons exhibit a larger V2 than the mesons. While 
the pion and kaon V2 reach a similar maximum of V2 ~ 0.14 at pT ~ 2.5 GeV/c, 
the baryon V2 continues to rise until it reaches a maximum of V2 ~ 0.20 at px ~ 
4.0 GeV/c. For still larger pT, the V2 values exhibit a gradual dechne until V2 for all 
particles is consistent with V2 ~ 0.10 at pt ~ 7 GeV/c. Fine detail cannot yet be 
discerned at > 7 due to statistical and systematic uncertainties. At these higher 

IQ CI 

Pt values one expects that the dominant process giving rise to V2 is jet-quenching^ 
where hadron suppression is larger along the long axis of the overlap region than 
along the short axis"" o' oo. For very large energy loss, the value of V2 should be 
dominated by the geometry of the collision region. Fig. [16] shows a comp aris on of 
V2 datcP^ for 3 < pt < 6 GeV/c compared to several geometric modelJ^Sl. This 
comparison seems to indicate that V2 in this intermediate pt range is still too large 
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to be related exclusively to quenching. 




b (fm) 

Fig. 16. The four-particle cumulant «J2 at 3 < pt < 6 GeV/c versus impact parameter, b, 
compared to models of V2 based on geometry alone. Data in this region does not agree well with 
these surface emission pictures. 



In the higher px regions, significant azimuthal structure will arise from jets, so 
non-flow correlations are thought to be significant in this regiori^. These effects 
have been studied in several ways. The four-particle cumulant V2 has been studied as 
a function of pt and the ratio of the four- a nd two-particle cumulants W2{4} /w2{2} is 
found to decrease with increasing ptEHII]^ This decrease is identified with a gradual 
increase in the contribution of jets to U2{2} (see Fig.[T7]left panel). The four-particle 
cumulant suppresses contributions due to intra-jet correlations but the statistical 
errors of the measurement are larger. One can also suppress jet structure in the V2 
measurement by implementing a Ary cut in the pairs of particles being used in the 
analysiJ^. In this high Pt particle is correlated with other particles in the 

event that are separated by a minimum A77. This method relies on the assumption 
that jet correlations do not extend beyond a given Ary range. Interactions of jets 
with the medium in nuclear collisions however can change the structure of jets and 
extend the correlations in Ary beyond the widths observed mp+p collisions^. This 
method therefore is not guaranteed to eliminate non-flow from jets. The problem 
of measuring V2 without non-flow and of measuring modifications of jet structure 
by the medium are entirely coupled. If one is known, the other is trivial. 

Other methods for suppressing non-flow include measuring correlations between 
particles at mid-rapidity and and an event-plane determined from particles observed 
at forward rapiditjE31. In the extreme and the most effective case, the event-plane 
was reconstructed from spectator neutrons in a Zero-Degree Calorimeter to mea- 
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Fig. 17. The left panel shows the ratio t)2{4}/v2{2} for the 0%-70% centrality interval in 200 
GeV Au-|-Au collisions. The ratio falling below unity indicates the importance of non-flow and V2 
fluctuations. The larger reduction at high pT seems indicative of an increase in non-flow due to 
jets. The right panel shows the centrality dependence of the t)2{4} and V2 measured with respect 
to the event plane in 130 GeV Au-(-Au collisions. The inset shows the ratio of the two. 



sure V2 of produced particles near 77 = 0. An extension of analyses based on the 
change in correlations across various rapidity intervals is the analysis of the two 
dimensional correlation landscape for two-particle correlations e.g. (P N / A(j)Ai^^. 
After unfolding the two particle correlations one can attempt to identify various 
structures with known physics such as jets, resonance decay, or HBT based on their 
width in 77 and 0. The remaining cos(2A(/)) structure can then be used to estimate 
{^2)"^ + c^j- This method will be discussed below. 

As we study progressively smaller systems the connection between the nucleus- 
nucleus reaction-plane and the azimuthal structure breaks down. In the limit that 
one proton from each nucleus participates in the interaction, the reaction-plane 
defined by the colliding protons will not necessarily be related to the reaction-plane 
defined by the vector connecting the centers of the colliding nuclei. In order to 
facilitate a comparison between the px dependence of azimuthal correlations in 
large systems and small systems, the scalar product (uQ*) is used where u = e^'^'^ 
and Q* = e~'2*E2]. The mean of uQ* therefore yields a quantity that depends on 
{v2V2{pt)) and non-flow as follows: 

{uQ*) = (S, cos(7i(0p^ - $,))> = M{v2V2iPT)) + MS2, (5) 

where M is the multiplicity used in the sum. Fig. 1181 shows this quantity for p + p, 
d+Au, and three Au+Au centrality intervals. The p+p and c?-|-Au data are repeated 
in each of the panels. The most peripheral Au+Au collisions are shown in the left 
panel. The centrality bin shown is not usually presented since trigger inefficiencies 
for low multiplicity events makes it difficult to define the actual centrality range 
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Fig. 18. Second harmonic azimuthal correlations in p + p, d+Au, and Au+Au collisions. The 
quantity (S^ cos{n(0py — ^i))) = M{v2V2{pt)) + M52 facilitates comparisons between different 
systems particularly where a reaction-plane may not be well defined. 



sampled. In this case, the data has been published in order to compare uQ* between 
the most peripheral sample of events and p + p cohisions. The data in Au+Au has 
a similar shape and magnitude as the data in p + p. This suggests that peripheral 
collisions are dominated by the same azimuthal structure as p + p collisions; an 
observation consistent with two-particle Ary, A(f) correlation^^. The data from mid- 
central Au-|-Au collisions shown in the middle panel however, exhibit a magnitude 
and shape clearly different than p+p collisions. While uQ* for p+p, d+Au and very 
peripheral Au-|-Au collisions rises monotonically with pT, for mid central Au-|-Au 
collisions, the data rises to a maximum at pt = S GeV/c and then falls. For central 
collisions shown in the right panel, a similar feature is seen with data rising to 
a maximum at = 3 GeV/c, then falling until px = & GeV/c, where it begins 
rising again. This second rise is presumably a manifestation of non-flow at high 
Pt in central collisions. These data suggest that azimuthal structure in Au-|-Au 
collisions above pt = 6 GeV/c is dominated by jets. This is also consistent with the 
conclusions reached by examining the particle type dependence of V2 and Rc^^- 
In Fig. [18] the p + p data is replotted in each panel to facilitate a comparison 
between the shape and magnitude in p + p to that in Au+Au. In the absence of 
jet-quenching however, non-flow at high px is expected to scale with the number 
binary nucleon-nucleon collisions A'^binary The plotting format in Fig. [18] on the 



June 4, 2009 18:40 WSPC/INSTRUCTION FILE qgp4-v2chapter 



Review of elliptic flow V2 21 




Pj (GeV/c) 

Fig. 19. The left panel reproduces the data in Fig. 1181 The right panel shows R^^*^- Binary 

scaling of S2 from p + p collisions corresponds to R^Xa*^ = 1. In comparison to binary scaling, 
non-flow above pj- ~ 5 GeV/c in Au+Au collisions is suppressed by a factor of 5. Also shown for 
comparison is Raa for single particle spectra of charged hadrons. 



other hand, assumes that ^2 oc 1/M rather than Nbin/M'^ as would be expected 
for hard scattering. The multipU city h as been shown to scale as (1 — a;/iarrf)-^part + 
XhardNunary with Xhard ~ o_iil 97 | 98 | -pj^jg jg j-gferrcd to as the two-component 
model. In order to compare azimuthal structure in Au+Au collisions to A'binary 
scaling of p + p collisions we can form a ratio in analogy with Raa for single 
hadrons: 

j^uQ" ^ {uQ*)aa/Maa 

where Maa and Mpp are the multiplicities in A + yl and p + p collisions with AIaa 
taken according to the two component model. In the case that jet production in 
Au-|-Au collisions scales with the number of binary collisions, as hard processes 
are expected to, should be unity. The right panel of Fig. [19] shows R^aa ^'^^ 

charged hadrons in 0% — 5% central Au-|-Au collisions. For comparison, Raa from 
single particle charged hadron spectra is also shown in the figure, i?^*^ first rises 
abruptly with px to a maximum of 2 at px ~ 0.5 GeV/c and then falls to a value of 
0.25 at pt ~ 5 GeV/c. At px > 5 GeV/c i?^*^ is similar to Raa- This shows that 
jet-quenching suppresses the charged hadron spectra, and the azimuthal structure 
by a similar amount; confirming that the single hadron suppression is indeed related 
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to jet-quenching, i?^^ is complimentary to studies of Iaa^ the ratio of dihadron 
correlations in Au+Au and p + p collisionJ^. 

We note the presence of what appears to be a local minimum and local maximum 
Pt ~ 1-5 and 2.0 GeV/c respectively. It is not clear if this is a real feature or simply 
an artifact largely caused by the shape of the p + p data. In the case that it is a 
real feature, it is possibly related to the changing particle composition in Au+Au 
collisions where baryons with larger V2 values become more prominent. At pt = 
3 GeV/c, baryons and mesons in p + p collisions are created in the proportion 1:3 
while at the same pT in central Au+Au collisions the proportion is approximately 
1:1. i?^*^ will be an interesting quantity to investigate for identified particles. One 
can anticipate a quark number dependence at intermediate pt as seen in Rc p and 

V2- 



2.3. Multiply strange hadrons and heavy flavor 

The build-up of space-momentum correlations throughout the collision evolution 
is cumulative. Information about space-momentum correlations developed during a 
Quark-Gluon-Plasma phase can be masked by interactions during a later hadronic 
phase. For studying a QGP phase, it is useful to use a probe that is less sensitive 
to the hadronic phase. Multi-strange hadrons have hadronic cross-sections smaller 
than the equivalent non-strange hadrons, and the V2 values measured for hadrons 
such as (^-mesons (ss) and f2-baryons (sss) are therefore tho ugh t to be more sensi- 
tive to a quark-gluon-plasma phase than to a ha dronic phase^. 

Fig. [20l shows V2{pt) for the 0-mesorl22linn]_ 'pjjg rises with pT and reaches 
a maximum of approximately 15% at pt near 2 GeV/c. At intermediate pt, the (p- 
meson V2 appears to follow a trend similar to the other meson K'g. This observation 
suggests that either the 0-meson cross section is larger than anticipated or re- 
scattering during the hadronic phase does not contribute significantly to V2. The 
latter possibility requires that V2 is established prior to a hadronic phase, suggestive 
of development of V2 during a QGP phase. 

Measurements of multiply strange hadrons are interesting because they should 
be less coupled to the matter in a hadronic phase and therefore a better reflection of 
the QGP phase. Heavy quarks on the other ha nd {e.g. c harm and bottom quarks) 
may be less coupled to even the QGP matteJiSlMl. 

It's not a priori obvious 

that heavy quarks will couple significantly to the medium and be influenced by its 
apparent expansion. The extent to which they do couple to the medium should be 
reflected in how large V2 for heavy flavor hadrons becomes and how much the nuclear 
modiflcation {Raa) deviates from unity. Precision measurements of Heavy Flavor 
mesons or baryons are not yet available from the RHIC experiments. As a proxy 
for identifying D-r nesons, t he STAR and PHENIX experiments have measured non- 
photonic 

electronp™. 

Non-photonic electrons are generated from the weak- 
decays of heavy flavor hadrons and after various backgrounds have been accounted 
for can, with some caveat^^^, be used to infer the Raa and V2 of D-mesons. 



June 4, 2009 18:40 WSPC/INSTRUCTION FILE qgp4-v2chapter 



Review of elliptic flow V2 23 



0.25 
0.2 

0.15 
0.1 

0.05 




Au+Au 200 GeV 

□ A • (|) 

o K° 



2 3 

(GeV/c) 



4 5 



Fig. 20. V2(pt) for the </> meson. The (/> is composed of an ss pair and is expected to have a 
smaller hadronic cross-section than non-strange, or singly-strange hadrons. The (j> V2 is compared 
to Kg and A V2. 



The top panel (a) of Fig. [211 shows Raa for non-photonic electronJl^SIinZ! Prior 
to the measurement of non-photonic electron Raa, it was expected that heavy- flavor 
hadrons would be significantly less suppressed than light flavor hadrons. These ex- 
pectations based on a decrea se in the coupling of charm quarks to the medium 
because of the dead-cone effecPH, are contradicted by the data; At ~ 5 GeV/c, 
non-photonic electrons are as suppressed as pions. This suppression suggests a 
stronger than expected coupling of charm quarks to the medium. This coupling 
apparently also leads to significant V2 as seen in Fig. [21] (b). 

A lso sh own in the figure is a calculation of V2 and Raa based on a Langevin 
modePSI. In that model, the strength of the energy loss and momentum diffusion 
of charm quarks is characterized in terms of a diffusion coefficient (D). Raa and 
V2 for charm quarks is then computed for several values of D. Two of these values 
are shown in Fig. 1211 Although neither curve provides an entirely satisfactory si- 
multaneous description of V2 and Raa, the comparison suggests that the diffusion 
coefficient is large. This comparison only achieves rough agreement, but the calcu- 
lation illustrates the sensitivity of heavy flavor hadrons to transport coefficients of 
the QGP and V2 is an important quan tity t o measure for these hadrons. This is also 
in agreement with H. van Hees, et al^^^ where coalescence at the hadronization 
phase boundary is also considered and found to help improve the agreement with 
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Fig. 21. iJ^A and t)2 for non-photonic electrons. The two measurements can be used in conjunction 
to provide constraints on models. Data are compared to a model with the conjectured minimum 
viscosity-to-entropy ratio rj/s and one with rj/s four times that value. 



data. 

2.4. Fluctuations and Correlations 

Comparisons between data and models are complicated by uncertainties in the 

initial eccentricity and by uncertainties in the data. Estimating transport quantities 

from the data may require a precision comparison between eccentricity and V2 so 

it is important to reduce the uncertainties in both. As discussed in the previous 

sections, a CGC model of the initial conditions yields eccentricity values typically 

30% larger than a Gluaber model while the fluctuations ((Tg) are still of the same 

width . The ratio of a^/e in a CGC model is therefore smaller than in a Glauber 

One can expect the statistical fluctuations in eccentricity to show up as 

dynamical fluctuations in V2 measurements. Measuring the dynamic V2 fluctuations 

in conjun ction with (V2) ca n therefore provide an additional constraint on the initial 
conditioniiMlMMlal, 
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Several methods have been employed for measuring V2 and the vari- 
ous metho ds have different depende ncies on non-flow correlations and V2 
fluctuationjM1114|115|116|71|59|117|ll8l. 

The differences between these measure- 
ments give information on non-flow correlations and V2 fluctuations. If one uses a 
two particle correlation to estimate V2, then one finds V2{2}'^ = (cos(2(0i — 4>j))) — 
(^^2)^ + +S2 where the average is over all unique pairs of particles. V2 is the single 
particle anisotropy with respect to the reaction plane V2 — {cos{2{(f> — 5"))) and 62 
is the non-flow parameter which summarizes the contributions to (cos(2((/)i — (f)j))) 
from correlations not related to the reaction plane. If one uses a 4-particle cumu- 
lant W2{4} calculation, then for most cases the non-flow term will be suppressed 
by large combinatorial factors and V2 fluctuations will cont ribute with the opposite 
sign. For Gaussian fluctuations, f2{4}^ ~ {^2)'^ ~ '^v^^- Without knowing S or 
cr^2, one cannot determine the exact value of {v2)- Rather, (w2)^ could lie anywhere 
between W2{4}^ and {v2{2}^ + V2{'i}^)/2. 




| ^ ^ ; ^ ; ^ 

2 4 6 8 10 12 14 
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Fig. 22. The quantity a^^^ = 52+2a'^^ estimated from the from the difference tJ2{2}^ — 1J2{4}^. Also 
shown is a parametrization of V2{1D}^ — V2{2D}^; the difference between (cos(2A0)} extracted 
for all correlations and that extracted by first fitting other structures in the correlations identified 
with various non-fiow sources. 



It is advantageous to confront various models with the data that is experimen- 
tally accessible. The difference between the two- and four-particle cumulants in the 
case of Gaussian V2 fluctuations is: 

«2{2}2-w2{4}2«<52 + 2<. (7) 

The term S + 2a'^^ is also approximately equivalent to the non-statistical width 
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of the distribution of the length of the flow vector distribution {dN/d\q2\) and is 
called (TjQj. The flow vector for the n*-^ harmonic is defined as qn.x = Sicos(n0i) 
and qn^y = sm{n(j)i). Fig. [22l shows a^^^ extracted from the difference between the 
two- and four-particle cumulantjm. In the case of Gaussian fluctuations, higher cu- 
mulants such as W2{6} are equal to ti2{4}. In this case, the quantity (7^^^ and V2{2\^ 
summarizes the information available experimentally from the second harmonic 
flow vector distribution. No more information can be accessed without applying 
more differential techniques or by making assumptions about the shape or cen- 
trality dependence of flow, non-flow, or flow fluctuations. An example of a more 
differential analysis is al so sh own in Fig. [221 where two particle correlations have 
been fit in A<^-Ar7 spacJ^^^. Terms identified with various non-flow sources have 
been included with the fit and the remaining cos(2A0) modulation is then identified 
as V2{'2D}^. In the case that the sources of non-flow are correctly parametrized, 
S = V2{ID}^ — V2{2D}'^ , where V2{1D}^ is (cos(2A(/))) integrated over all azimuthal 
structure. Then f2{2_D}^ = {v2)^ + cr^^ - This procedure is discussed below. 



■J— nucleon participants 

■*- constituent quark participants 

■-■ color glass (fKLN) 




impact parameter (b) (fm) 



Fig. 23. An upper limit on the ratio cr^^/{v2) is shown. The upper Umit can be established from 
(Ttot and t)2{4}. The upper limit is compared to various models of eccentricity fluctuations. The 
existence of non-flow correlations implies that the true value of cr^j / {v2 ) may be significantly below 
the upper limit, potentially challenging the Montc-Carlo Glauber calculation of the eccentricity. 



Even without attempting to disentangle flow fluctuations from non-flow cor- 
relations, the assumption that non-flow is a positive quantity (consistent with 
V2{ID}^ — V2{2D}^) can be used with (7^^^. to provide an upper limit on V2 fluctu- 
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< < -f- (8) 

To facilitate a comparison between this limit and models of the initial eccentricity 
the upper limit on is compared to jf^ in the model. To form this ratio appro- 
priately, the same assumptions should be made for a^^ ^^id {V2) i.e. zero non-flow. 

< ^ ^ V2{2V-V2{'iV 

(«2>2 («2{2}2 + «2{4}2)/2 V2{2}^ + V2HV ^' 

This upper limit is shown in Fig. [23] and compar ed to several models of ag/e. The 
models include two Glauber Monte Carlo one using the coordinates of 

participating nucleons to calculate the eccentricity, the other using the coordinates 
of constituent quarks confined insid e the nucleons. The constituent quark Monte 
Carlo Glauber Model (cqMCGj^l 

treats the nucleus as 3 x A constituent quarks 
grouped in clusters of three confined to the size of a hadron. This increases the 
number of participants by roughly a factor of three, reducing the fluctuations in ec- 
centricity. The correlations between the constituent quarks required by confinement 
partially counteract this effect since those correlations act to broaden the eccentric- 
ity distribution. The net effect, however, is a narrowing of the distribution. Also 
shown is a Color Glass Condensate (CGC) based model which yields eccentricity 
values 30% larger than the Glauber models leading to a reduction of (Je/e. The 
Monte Carlo Glauber model based on the eccentricity of nucleons already exhausts 
most of the width a^^^ — 62 + Zcr^^. This shows that the statistical width of the 
eccentricity fluctuations in the Glauber model already accounts for almost all of 
the non-statistical width of the flow vector distribution thus leaving little room for 
other sources of fluctuations and correlations. This is particularly challenging since 
non-flow has been neglected in setting the upper limit and the only V2 fluctuations 
considered are those arising from eccentricity fluctuations. We have therefore ne- 
glected fluctuations that would arise during the expansion phas e ■'. One 
can write the total width including these terms: 

e 

where the middle term in the right-hand-side is the V2 fluctuations from eccentric- 
ity fluctuations and the final term is V2 fluctuations from the expansion phase. The 
middle term arises from the approximation that to first order cry^/v2 = o^je. This 
approximation is prevalent in the literature. The last term can be related to the 
Knudsen number of the matter during the expansior l^^^ i Measurements demon- 
strating the existence of non-flow or dynamic V2 fluctuations {(7^2. dyn) therefore 
would challenge the model. The CGC or cqMCG models provide a more likely de- 
scription since a^/e is smaller in those models than the upper limit on a^^lk^'i)- 
The upper limit on (Jt,^ provided by provides a valuable test, therefore, for 
models of the initial eccentricity and can help to reduce the uncertainty on e; an 
essential component in extracting meaning from the value of V2. 



= ^2 + 2(«2-)' + 2<,,,„, (10) 
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2.4.1. Two Dimensional Correlations and V2{'iD} 




2Nbin/N part 2Nbin/N part 



Fig. 24. The top panels show simulated data for the two-particle correlation measure Ap/ ^Pref- 
The simulations are tuned to yield data similar to p + p collisions (left) and Au+Au collisions 
(right). The bottom left panel shows V2 estimated from the 2- and 4-particle cumulants. The 
shaded regions show the V2 derived by integrating over all structure in Ap/^py-e{ (upper shaded 
region) or by performing a 2-Dimensional fit to the surface of Ap/^/p^ (lower shaded region). The 
bottom right panel shows the 2- and 4-particle cumulant data transformed to compare directly to 
Ap[2]/7ft;?: {2} = nv2{2}^ and {4} = nvii'l}'^ . 

One way to study non-flow contributions to two -particle correlations is to mea- 
sure the correlations as a function of Arj and A<^^. This allows different sources of 
two particle correlations to be studied where each source is identified by its char- 
acteristic dependence on A(f) and A77. Additional information can be obtained by 
including information about the charge-sign dependence of the correlations. An ex- 
ample of such an analysis is shown in Fig. 1241 The figure displays four panels. The 
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top panels show the correlation density 

^ = ^^, (11) 

y/Pief \/Picf 

where p is the pair density and prof is the product of the single particle densities. 
This normalization is chosen to search for deviations of the correlations in large sys- 
tems from those in small systems. If Au+Au collision were simply a superposition 
oi p + p collisions for example, Ap/^pi-cf would be the same in p + p and Au+Au 
collisions. This figure was produced based on simulated data which was tuned to 
match real data in very peripheral A u+Au or p+jP colhsions (top left) and 20%-30% 
central Au+Au collisions (top rig The correlations in these two systems 

have been found to be very different. The p + p collisions exhibit structures charac- 
teristic of fragments from string breaking (a narrow ridge at Arj — independent of 
A(f)) and fragments from semi-hard scattered partons or mini-jets. These mini-jets 
yield a two-dimensional Gaussian correlation at 0, 0, and a broad ridge at Acj) = tt. 
The away-side jet can sweep over a wide Arj range since the partons can have a 
momentum within the proton or Au nucleus. For semi-central and central Au+Au 
collisions, the correlation landscape is drastically different with the most prominent 
feature being V2 giving rise to a clear cos(2A0) shape. If the shape of the various 
non-flow terms structures is well understood, the correlation landscape can be fit 
and (v|) can be extracted, independent of the sources of non-flow. This procedure 
depends on having an accurate description of the shape of the non-flow sources. 
Some of these are easily identified based on their charge dependence or their char- 
acteristic shapes. Other sources may be less easily identifiable though, particularly 
if they become modified by the medium in Au+Au collisions. 

The bottom panels of Fig. [24] show a proof-of- principle extraction of (ul) based 
on the simulated data. The left panel shows while the right panel shows the 

per-particle measure Ap[2]/ ^fp^. Also shown are the two-particle and four-particle 
cumulant data. The upper hatched region in the bottom left panel shows 

{vD = ^Ap[2]/^ (12) 
n 

extracted without separating out non-fiow; called V2{1D}. n is the multiplicity 
of measured tracks. The lower hatched region shows the same without including 
the non-flow structures in the calculation; called d|{2D}. Data are plotted versus 
2-/Vbinary/A^part SO that most Central collisions are on the right. This procedure makes 
use of the two-particle correlations landscape to separate different contributions 
to the azimuthal structure. The fit procedure does require some assumptions be 
made in order to separate non-fiow from v These include that V2 and (T ^2 are 
independent of A77 and that higher harmonics of w„ and a^^ do not contribute. 
Relaxing those assumptions may make it difficult to distinguish between non-flow 
and more complicated correlations related to the reaction plane. More information 
can be made use of however, by examining the charge, sign dependence and particle- 
type dependence of the various correlation structures. This is therefore a promising 
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method for disentangling V2 and non-flow. 



2.5. Scaling Observations 

Elliptic flow measurements represent an extensive data set. V2 has been measured 
for 0.1 < pt < 12 GeV/c, for —5 < 77 < 5, for mesons from the pion to the (j) and 
J/V', for baryons from the proton to the 51, and for transverse particle densities 3 < 
< 30. And yet, given the complexity of heavy-ion collisions and such a large 
data-set, the measurements exhibit many surprisingly simple features. These we can 
summarize in terms of simple scaling observations where a large amount of data is 
found to behave in a regular and simple way when plotted versus the appropriate 
variable. The observation of a particular scaling then motivates the question: why 
does the data only depend on that variable? These scaling observations, therefore, 
not only allow us to summarize large amounts of data in a simple form, but they also 
suggest simple physical explanations for the data with perhaps deeper implications. 
In this section I review several observed scaling laws. 



2.5.1. Longitudinal scaling 

Fig. [T2I shows the centrality dependence of V2{'q)^^ ■ Although more detailed mea- 
surements at small 77 show that V2 is approximately independent of 77 for |?7| < ff^, 
the data extending to larger rj exhibit a nearly triangular shape: having a maximum 
at 77 = with a nearly linear decrease with I77I. A similar shape is seen for ^s„„ = 
200, 130, 62.4, and 19.6 GeV Au-hAu colhsionJI^l. 

Fig. [25l shows V2{\'q\ —yteam) for ^s„„ from 19.6 to 200 GeV; one order of mag- 
nitude in ^s„„. The data are for the 40% most central collisions. Ideally the x-axis 
would display y — ybeam but data on identified particle V2 spanning such a large 
range of rapidity are not available. One finds that within errors, all data lie on a sin- 
gle curve. This suggests a smooth variation of the development of space-momentum 
correlations from forward rapidity to mid-rapidity. This scaling observation also im- 
plies that the value of V2 obtained at mid-rapidity is a smooth function of yteam 
or equivalently of log(^s„„); consistent with the smooth trend seen in Fig. [3] for 
V2 above ^s„„ of approximately 5-10 GeV. An energy scan at RHIC extending 
down to v^s^viv — 5 GeV will make it possible to investigate this trend with better 
precision and with a single detector, eliminating many systematic uncertainties 
This simple trend may be confirmed with more precision or perhaps deviations will 
point to a softest point in the equation-of-statJI^Sn^. 



2.5.2. Kinetic Energy and Constituent Quark Number Scaling 

At low pt, V2 is or dered b y mass with heavier particles having a smaller V2 value 
at a given 

valuOT 

This ordering is indicative of particle emission from a 
boosted source with the boost larger in the in-plane direction than the out-of-plane 
direction. Indeed, blast-wave fits implementing this scenario agree very well with 
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Fig. 25. V2{'ri — ybeam) for ^ variety of energies. The centrality interval used is 0% — 40%. Plotted 
in this format, data at all energies fall on one curve. This scaling is also seen in dN/dy{r) — ybeam) 
and is referred to as limiting fragmentation. V2 is largest at middle rapidity (»; = or r; — yteam = 
-Vbeam) and vanishes at = ybeam iv - ybeam = 0). 



the data in this region^. It is also found that in thi s sam e pt region, when V2 is 
plotted versus jtlt — rriQ all data fall on a common rriT — Too is the particles 

transverse kinetic energy and sometimes labeled KEjC^ 

Fig. [26] shows V2 versus tot — toq for particles ranging in mass from the pion 
with mass of 0.1396 GeV/c^ to the S with mass of 1.321 GeV/c^. The measurement 
is made for the 0-80% centrality interval in 200 GeV Au+Au collisions. Similar 
scaling has also been demonstrated for 62.4 GeV collisionJ^ll. The data exhibit 
obvious trends. At low to^ — toq, V2 values for all particles rise linearly with no 
apparent differences between the particles with different masses. Near tot — toq = 
0.8 GeV/c^, f2(TOT — Too) for mesons and baryons diverges. The meson V2 begins to 
saturate, obtaining a maximum value of 14-15% near uit — toq — 2.5 GeV/c^. The 
baryon V2 continues to rise, obtaining a maximum value of approximately 19-20% 
at Tot — Too = 3 to 3.5 GeV/c^. The relative masses of the baryons and mesons 
do not seem to be relevant, rather the number of constituent quarks in the hadron 
determines the V2 values in this range. The mass dependence can be better checked 
using the 0-meson which has a mass slightly larger than that of the proton. The 
statistical significance of the cj) V2 is limited but measurements seem to indicate that 
the 4> lies closer to the mesons than to the baryons i.e. closer to the particles with 



a common number of constituent quarks than to particles with a common mass^. 
The observation of the quark-number dependence of V2 at intermediate px led 
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Fig. 26. V2 for a variety of particles plotted versus my — mo where m|, = + rrig, and mo is 
the rest mass of the particle, mj- — mo is also the transverse kinetic energy of the particle KEj-. 



to speculation that hadron formation through the coalescence of dressed quarks 
at the hadronization phase boundary could lead to an amplification of V2 with 
bar yons Retting amplified by a factor of 3 while me sons were amplified by a factor 
of jl30|131|132ll40|133|134|135|136|137|138|139|140| ^his picture was subsequently 

stre ngthene d by the observation that a similar quark-number dependence arises in 
the ratio of the single particle spectra in central collisions to that in 
peripheral collisions. At intermediate pr the Rcp values for various particle species 
are also grouped by the number of constituent quarks, with baryons having a larger 
Rcp- The larger Rcp for baryons signifies that baryon production increases with 
collision centrality faster than meson production; an observation consistent with 
the speculation that hadrons from Au+Au collisions are formed by coalescence 
such that baryon production becomes easier as the density of the system increases. 
The more general and less model dependent statement is that the baryon versus 
meson dependence arises from high density and therefore most likely from multi- 
quark or gluon effects or sometimes called" higher twist" effects. The combination 
of large baryon V2 and large baryon Rcp also immediately eliminates a class of 
explanations attempting to describe one or the other observation: e.g. originally it 
was speculated that the larger Rcp for baryons might be related to a smaller jet- 
quenching for jets that fragment to baryons than for jets that fragment to mesons. 
This explanation would lead to a smaller baryon V2 and is therefore ruled out by 
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the larger V2 for baryons. The same can be said for cofor transparency 

which would account for the larger baryon Rcp in this px region but would predict 

a smaller baryon V2. Color transparency may still be relevant to the particle type 

dependencies at px > 5 where Rcp for protons is slightly larger than Rcp for 
1 4^ 

pions and the V2 measurements are not yet precise enough to conclude whether 
the baryon V2 is also smaller than the meson V2- This is a topic that needs to be 
studied further. 



Vj/n 
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A Kg(n=2) • A + A(n=3) 



A* 



12 3 

Transverse Momentum p^/n (GeV/c) 



Fig. 27. V2 for Kg and A scaled by the number of their constituent quarks (n) and plotted versus 
px/n. The data appear to fall on a universal curve which has been taken as an indication of 
hadron formation via coalescence of quarks from a flowing medium. 



In a coalescence picture, the final momentum of the observed hadron would de- 
pend on the momentum of the coalescing constituent quarks. The exact dependence 
is not known but a relatively good scaling of V2 for Kg and A was found when V2 jn 
was plotted as a function of pt /n. Such a scaling implies that the momentum of the 
hadron is simply the sum of the momenta of the coalescing quarks. Fig. [27] shows 
V2/n versus px/n for X^-mesons and A-baryons. The scaling appears to be good 
throughout the whole pT range but part of this perception is due to the decrease 
of V2 for both particles at small px- When a ratio is taken between the V2 / n{pT / n) 
values, a clear deviation from scaling is seen in the lower pT region. A combination 
of the Tot — Too scaling in Fig. [26] and the V2/n scaling in Fig. [27] will lead to a 
good scaling over the whole measured momentum range; since U2(tot — toq) for all 
particles fall on a single line at low — toq , dividing the x- and y-axis by n will 
not destroy that scaling seen in Fig. [26] Plotting V2/n versus (tot — mo)/n should 
therefore provide a good scaling across a large kinematic range. 

Fig. [28] shows V2/'n versus (tot — niQ)/n for 200 GeV Au+Au collisions in three 
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Fig. 28. A more detailed study of quark number scaling in Au+Au collisions. In the left panels (a, 
c, and e) v^jn is shown versus (my — mo)ln for three centrality classes. Hydrodynamic models 
are also shown for comparison. Data are fit to a single curve. In the right panels (b, d and f) the 
ratio of the data and hydro model to the fit function are shown. 



different centrality interval^SQI Data for i^Tg-mesons, A-baryons and S-baryons are 
sliown. The left panels show the data with a hydrodynamic calculation and a fitting 
function. The phenomenologically motivated function 

a + fox + cx^ a 
1 + exp[ ' ] ^ 

with X — (niT ~ mQ)/n, describes the data well for the three centralities. The 
function captures the rise then saturation and steady decline seen in the data. We 
ascribe no physical meaning to the function or the five fit parameters but simply 
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use it as a convenient reference. The right panels show the ratio of the data and 
the hydro model to the fit function. For reference, the fit parameters are shown in 
Table [TJ The data is in good agreement with the fit function for all centralities while 
this hydro model calculation does not agree well with the data in any centrality. 







Fit parameters for Eq. [T3l 




Centrality 


a 


bed 


e 


40%-80% 


20.0e-02 


0.0 0.0 -1.19e-02 


2.37e-01 


10%-40% 


16.4e-02 


-4.53e-03 0.0 2.61e-02 


2.40e-01 


0%-10% 


8.96e-02 


-4.08e-03 0.0 6.52e-02 


2.70e-01 



Table 1. The fit parameters describing tlic curves in Fig. 1281 



There is a systematic deviation from the ideal n scaling at (tot — mo)/n > 
0.8 GeV/c^ with Kg mesons having slightly larger V2/n values than A baryons. 
This deviation from ideal scaling was predicted based on the inclusion of higher 
fock st ates in th e hadrons or the inclusion of a finite width in the hadron wave 
functioiPlSSl. Deviations can also arise in a hadronic phase when the hadronic 
cross sections are relevant. In the case that hadronic cross-sections are an impor- 
tant factor, higher statistics data for baryons and (j> mesons should deviate from 
their respective groups. We also note that hadronic cascade models also obtain ap- 
proximate V2 /n s caling due to the use of the additive quark model for hadronic 
cross- section a ^ '^^[ On the other hand, these models under-predict the integrated V2 
by a factor of two. We also note that non-fiow contributions can affect the scaling 
observed in this range and the particle-type dependence of non-flow sources is still 
being investigated. 

In Fig. [29] we investigate the breaking of ideal scaling in more detail with data 
integrated over a larger centrality interval. While this reduces the statistical un- 
certainty, it also introduces uncertainties due to the large centrality bin width. In 
particular, when particle yields have different centrality dependencies, the average 
eccentricity of events producing a particle can deviate from particle to particle. 
For example, the enhancement of baryons in central collisions will mean that the 
average baryon comes from a more central event than the average meson. Given the 
decrease of V2 with centrality, this can lead to a decrease of baryon V2 simply due 
to the wide centrality bin. Although there are caveats and systematic errors still to 
be quantified, we note that the baryons in Fig. [29] appear to lie systematically and 
significantly below the mesons. Self-similar curves are fit to mesons and baryons. 
The curves appear to describe the data. We note that the two self-similar curves 
shown in Fig. [29] can be nearly unified if we replace n with n+1. This demonstrates 
that the naive constituent quark scaling is violated to the extent that baryon V2 is 
actually closer to 4/3 the meson V2 rather than 3/2. The connection of the baryon 
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Fig. 29. V2/n versus (niT — mo) / n for minimum bias Au+Au collisions. The quark number scaling 
appears to be violated when integrating over a wide centrality bin. Mesons and Baryons are fit 
with the same functional form but with different parameters. The scaling is violated to the extent 
that n + 1 would give a better agreement i. e. baryon V2 is closer to four-thirds the meson V2 rather 
than three-halves. 



versus meson dependence and the number of constituent quark scaling appears to 
not be as directly connected to the number of constituent quarks as originally con- 
ceived. Whether this is indicative of higher fock states, the wave-function of the 
hadrons, an as yet un-accounted for experimental systematic error, or something 
else is yet to be determined. The systematic uncertainties based on the particle-type 
dependence of non-flow are still being investigated. 



2.5.3. System-Size Scaling 

The system-size dependence of V2 can be studied by looking at the centrality de- 
pendence of V2 or by colliding smaller nuclei. Ideal hydro predictions, having a zero 
mean-free-path assumption, should be independent of the system-size. In this case, 
given the same eccentricity, the V2 should be independent of system size. One can 
try to account for the change in eccentricity by dividing V2 by eccentricity from a 
model but this introduces a large amount of uncertainty. Another approach is to 
study the shape of V2{pt) to see if that varie The left panel of Fig. [30l shows 
V2 measured in Au-t-Au and Cu-|-Cu collisions for several centrality intervals In 
the right panel, V2{pt) is scaled by 3.1 times the mean V2 for that data set. 3.1(t;2) 
was taken as a proxy for the eccentricity of the collision system, and this proxy 
is not inconsistent with models of eccentricity which are quite uncertain. What is 
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Fig. 30. Left panel: V2 versus PT in Au+Au and Cu+Cu collisions for four centrality intervals. 
Right panel: the same V2 scaled by 3.1 times the integrated V2. The scaling demonstrates that 
the shape of 1)2 (pt) is approximately independent of centrality and system size. A function is fit 
to the Au+Au data and shown as a solid orange line. 



best demonstrated by this scaling, is that although the magnitude of V2 changes 
significantly for the different centralities and systems, the shape of V2{pt) is very 
similar. 

The invariance of V2{pt) with system-size can be taken as an indication that 
the viscosity of the expanding medium created in heavy-ion collisions can not be 
large when V2 is established; Large viscous effects should introduce a system-size 
dependence to ^a fPT'') with viscosity causing V2 to saturate at lower pT values in 
the smaller systeirU^^. Hyd rodynamic calculations including viscosity confirm this 
ide al •^^'^^ l ^'^^ l ^^^ l ^^^ l ^^^ l ^ ^. To look more carefully for a system-size dependence in 
the shape of V2 we plot the ratio of the scaled data to a curve fit to the Au-|-Au data. 
The results are shown in Fig. [3T1 The Cu+Cu data systematically deviate from the 
Au-I-Au data. The px dependence of the ratio indicates that the Cu-|-Cu data begins 
to saturate before the Au-|-Au data. This leads to a ratio that first rises then falls. 
This would happen the other way around if the Au-|-Au data saturated first. The 
uncertainties in the figure are large but the shapes are still significantly different. 
The system-size dependence of V2{pt) may be a valuable tool for estimating the 
viscosity of the matter created in heavy-ion collisions. 

Although the data on V2 includes many particle types, a wide kinematic range 
in pt and ry, a variety of system-sizes and a wide range in center-of-mass energy, 
we've been able to identify several regular features of the data. These include a 
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Fig. 31. The 'U2{PT)/3.1(f2) data from the right panel of Fig. 1301 scaled by a function fit to the 
Au+Au data. This figure illustrates that their seems to be a significant difference between the 
shape of 1)2 (pt) in Au+Au and Cu+Cu collisions with the Cu+Cu collisions exhibiting a more 
abrupt turn-over — i.e. Cu+Cu data first rises faster with pT then falls faster at pT > 1 GeV/c. 



nearly linear rise of V2 at mid-rapidity with log(y^J^): 

V2 = 0.008 + 0.0084 log(V^) (14) 

for 0%-20% central Au-f Au or Pb+Pb collisions and a pr, mass and particle- type 
dependence that can be parametrized by 

a + bx + cx^ a 

1 + exp[ ^\ ' \ ^ 

where x = [rnT — m^^jn^ while Viij]) to good approximation decreases linearly from 
it's maximum at mid rapidity to beam rapidity. This linear rise may be a trivial 
consequence of the log(^s„^) dependence of mid-rapidity vi or vice-versa. 

3. Confronting the Hydro dynamic Paradigm with RHIC Data 

We have discussed the hydrodynamic model extensively in this review as a con- 
venient reference for how well the matter produced in heavy-ion collision converts 
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spatial deformation into momentum space anisotropy. Hydrodynamic models of 
heavy-ion collisions have many uncertainties. These include, uncertain initial con- 
ditions, uncertain thermalization times, and uncertain freeze-out conditions. A suc- 
cessful description of data using a hydrodynamic model offers the promise of not 
only establishing the attainment of local equilibrium but also the promise of provid- 
ing information on the Equation-of-State of the matter and its transport properties. 
The uncertainty in the models, however, are large and it has not yet been possible to 
extract this desired information with satisfactory cer tainty. In addition, the possi- 
bility that significant V2 arises from initial-state effect could call into question 
the applicability of hydrodynamics and the need for prolific final-state rescattering. 
Measurements of tw o partic le correlations, which have often been interpreted as 
arising from mini- need to be reconciled with the idea of a locally ther- 

malized matter with extensive final-state rescattering. If the hydrodynamic models 
and data are irreconcilable, the paradigm will, of course, have to be abandoned. 




Fig. 32. The four panels show px spectra and V2 (top and bottom) for pions and protons (left and 
right). Data from 200 GeV Au+Au collisions is compared to a variety of hydrodynamic models. 
Most models do not agree with pT spectra and V2 simultaneously. 



To check for consistency with hydrodynamic model j^ ^ l '^^ ^ ^^^ l ^^^ l ^^^ l, the 
PHENIX collaboration created a comprehensive comparison between heavy-ion 
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data on pT spectra and V2/i^. The inclusion of a comparison to HBT data was 
hampered by the lack of predictions from some of the models. The comparison to 
Pt spectra and V2 is shown in Fig. 1321 The left panels show px spectra with pi- 
ous in the top panel and protons in the bottom. The right panels show V2 for the 
same particles. The combination of data on V2 and spectra provide a stringent test 
for the models as some models can reproduce one quantity but only by adjusting 
parameters in such a way that the agreement with other observables is spoiled. 

The models shown in the figure differ in several ways. Models that include a 
phase transition and a QGP phase are shown with solid lines while models without 
a pure QGP phase are shown as dotted lines. Including this phase transition acts to 
reduce the value of V2 since the equation of state is soft during the transition. This 
means that the speed of sound drops (in these models to zero), so that conversion 
of coordinate space eccentricity to momentum space anisotropy is halted during 
the phase transition. In the case that the models, do approximately match the pion 
spectra and V2: the most directly observable consequence of the lack of a phase 
transition is on the proton spectra and proton V2 ■ The proton spectra end up being 
too soft, and the splitting between proton and pion V2 is reduced with the proton V2 
becoming larger. This is somewhat counter-intuitive but is a consequence of fixing 
the parameters to match central data. 

The models also differ in their treatments of the final hadronic stage. The cal- 
culations from Teaney et al. include a hybrid model that uses a hadronic cascade 
(RQMD) for the final hadronic evolution. Hirano and Kolb do not use such an 
afterburner but allow the particle abundances to stop changing at a temperature 
above the temperature at which they stop interacting; chemical freeze-out happens 
before kinetic freeze-out. Huovinen on the other hand, maintains chemical and ki- 
netic equilibrium throughout the expansion. These different treatments have very 
important consequences for the particle-type dependence of the pr spectra and 
V2- Huovinen's treatment can reproduce the V2 for pions and protons, but only 
at the expense of under-predicting the number of protons; a direct consequence 
of maintaining chemical equilibrium until the final freeze-out at a relatively low 
temperature. 

The only model which compares well to all the data is Teaney's model including 
a QGP phase, a phase transition, and a hadronic phase modeled with RQMD. 
Such a hybrid model adds significantly to the number of tunable parameters as 
compared for example to Huovinen's model. On the other hand, the Teaney model 
shows that some particle types are less affected by the hadronic phase and therefore 
less sensitive to some of the uncertainty in freeze-out prescription. Fig. [33] shows 
the Teaney calculation with Hydro only versus Hydro-|-RQMD. The particle species 
least affected by the inclusion of a hadronic afterburner, are the 0-meson and the 
ri-baryon. This arises presumably from the small hadronic cross-section for these 
hadrons. This suggests high-statistics measurements for these particles are a viable 
way to avoid uncertainties in the effects of hadronic re-scattering. 
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Fig. 33. The inverse slope parameter Tgi^pe for a variety of particle species when a heavy-ion 
collision is modeled with only a hydrodynamic model compared to a hybrid model which includes 
a stage modeled with a hadronic cascade model RQMD. 



Besides the uncertainty in the freeze-out prescription, there is uncertainty on the 
eccentricity of the expanding fire-ball at the start of the conjectured hydrodyna mic 
evolution. Fig. [34] shows a hybrid hydro+cascade model compared to V2 

datJTH. 

Two model curves are shown: one with a Color-Glass-Condensate (CGC) initial 
eccentricitjESI, the other with a Monte-Carlo-Glauber (MCG) eccentricity. As dis- 
cussed previously the CGC eccentricity is larger than the MCG eccentricity; this 
leads to an over- prediction for V2- On the other hand, this hybrid model does not 
include viscous effects in the QGP phase so the difference between the hybrid-|-CGC 
prediction could be related to viscosity. In fact, since viscosity acts to reduce V2, 
the hybrid-|-MCG curve shows that there is no room for viscosity in this model. 
This viol ates t he lower bound on viscosity de rived based on quantum mechanical 
argument and also later from string theorjH^. Clearly, to estimate the viscos- 
ity allowed, or required by the data, the uncertainty on the initial conditions must 
be reduced. As discussed previously, the measured quantity a^^^ — 62 + 2cr,y^ pro- 
vides a sensitive test of the models of the initial conditions and needs to be carefully 
compared to the hydrodynamic model predictions with various initial conditions. 
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Fig. 34. A hybrid hydrodynamic model showing the uncertainty in the model calculations arising 
from two different models for the initial conditions. 



3.0.4. Transport Model Fits 

An approach to circumvent the uncertainties in the hydrodynamic models has 
been outhned in Ref.IISnilSIl where wj/e is fit as a function of The fit function 

is used to infer how close the data come to a saturated value in the collisions with 
the highest density achieved. The fit function is constrained by how V2/£ should 
approach the high density and the low density limits. One can construct different 
equations but in a transport code, th e foll owing was found to represent the approach 
to the zero mean-free-path limit welP^^l 

^ = ^ \ (16) 

where K is the Knudsen number and Kq is a constant of order one. Fig.[2S]shows the 
data and fit in the left panel and the inferred Knudsen number in the right panel. 
Based on this procedure it is found that RHIC V2 data are still some 20% below the 
saturation value anticipated within the fit function. This conclusion however, not 
only depends on the assumptions built into the transport model approach but also 
the centrality dependence of the eccentricity. The Color Glass Condensate model 
for example predicts a stronger centrality dependence for the eccentricity than the 
Monte Carlo Glauber model. As a consequence, this fit implies that if the initial 
conditions at RHIC are described by the CGC model, then the V2 data is closer to 
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Fig. 35. Left panel: V2 in Au+Au and Cu+Cu collisions scaled by eccentricity calculated in a 
CGC framework and plotted versus (1 / S){dN / dy) . The fit function and hydrodynamic limit are 
explained in the text. Right panel: The effective Knudsen extracted from the data and fit in the 
left panel. 



its saturation limit than if the MCG gives the correct description. This is counter 
intuitive and opposite to the conclusions reached based on real hydrodynamic cal- 
culations, which indicate that the lar ger CGC initial eccentricity allows more room 
for viscous effects in the QGP phase^SIISni xhe transport based fit circumvents 
the actual solving of hydrodynamics but the conclusions are dependent on the cen- 
trality dependence of the initial eccentricity which is strongly model dependent. 
The fit also includes the speed of sound as a free parameter. This effectively leads 
to an equation of state which has no phase transition but which is allowed to vary 
in the fit. A complimentary and perhaps better method for accessing the Knud- 
sen number and the viscosity is to study the shape change olmipT) for different 
system-sizes which avoids the uncertainty in the eccentricit y [ This is a work 
currently in progress. 



3.0.5. Viscous Hydrodynamics 

The apparent success of ideal hydrodynamic models to describe the gross features 
of RHIC data has led to the inference of small viscosity and the claim of the dis- 
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Fig. 36. Hydrodynamic calculations including viscous effects. The top panel shows 112 (pt) for the 
case that the initial conditions are described with a Monte Carlo Glauber model. The bottom 
panel is based on Color Glass Condensate initial conditions. The curves show results for different 
values of rj/s, the ratio of shear viscosity to entropy. 



covery of the perfect liquid at RHIC. The perfect liquid announcement was listed 
as the top physics story of 2005 by the American Institute of Physics and was 
widely covered in the popular press. Much recent work has gone towards includ- 
ing viscous effects in hydrodynamic calculati ons so that the viscosity can be more 
accurately estimateciMailM l50 | 151 | 152 | 153 | 

Fig. [36] shows one such calculation. 
The top panel shows results when the hydrodynamic evolution starts from Glauber 
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initial conditions while the bottom panel shows the case of CGC initial conditions. 
The results that come closest to the STAR datePIi are given by rj/s = 0.08 for 
the Glauber initial conditions and rj/s = 0.16 for the CGC initial conditions. The 
STAR non-flow corrected data are from Figure 4 of Ref.Sni. The 10%-40% central 
data from that reference was scaled to account for the difference between the 10%- 
40% centrality interval and the 0%-80% (minbias) centrality interval. The larger 
rj/s inferred based on the CGC model arises from the larger initial eccentricity 
which leaves more room for viscous effects that tend to reduce the V2- This contra- 
dicts the conclusions drawn from the transport model inspired fit, which allows the 
equation-of-state to change for the two different initial conditions. The pr depen- 
dence of the data is also better captured in the larger viscosity CGC scenario. The 
larger viscosity inferred from the CGC initial conditions gives a more pronounced 
turn over of V2{pt) which better describes the px dependence of V2. The comparison 
shown in Fig. [36] shows that hydrodynamic models including viscosity have a good 
chance of reproducing RHIC data as long as the shear viscosity to entropy ratio 
rj/s is less than ^ where is the conjectured lower hmit. 

3.0.6. Fluctuating Initial Conditions 

The comparison of V2 and other RHIC data to hydrodynamic models seems to in- 
dicate that when viscous corrections are included, a successful description of the 
data may be possible. There is uncertainty in this comparison, however, related 
to uncertainties in the initial conditions and in the freeze-out prescription. The 
uncertainty in the initial conditions can be addressed experimentally with mea- 
surements of V2 fluctuations which in turn require an understanding of non-flow 
correlations; The experimentally accessible information appears to reduce to V2{2}^ 
and V2{2}^ — V2{4:}'^ = S2 + 2a'^^. An alternative approach may be for hydrodynamic 
models to predict V2{2} and i'2{4} by including correlations and fluctuations in the 
models. Progress has been made in this direction. Early work relating to the effect 
of fluctuations in the initial conditions on hydrod ynamic calculations was carried 



fluctuations were indeed found to lead to V2 fluctuations as shown in Fig. 1371 Later 
it was suggested that correlations in the initial conditions could lead to w„ fluctua- 
tions of even and odd orders of n that would manifest themselves as non-s inusoida l, 
apparently non-trivial, two-particle correlations as seen in the RHIC datcPHH. 

Subsequent work following through on this idea shows that hydrodynamic mod- 
els with fluctuating initial conditions d o lea d to two-particle correlations with struc- 
ture beyond a simple cos (2A0) shapePS. The correlation structure arising from 
the fluctuations in the initial conditions is shown in the right panel of Fig. 1371 The 
model exhibits many of the features seen in the data including a jet-like peak, a 
near-side ridge, and an away-side ridge shifted away from = tt. All this structure 
arises without the explicit inclusion of jets in the model. The apparently exotic cor- 
relations do not appear in the model when a smooth initial condition is used. This 



out using the NeXSPheRIO 




eccentricity 
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calculation illustrates the importance of accounting for fluctuations in the initial 
conditions when interpreting the correlation landscape. It also demonstrates that 
complex interactions between jets and the medium, including mach-cones, are not 
needed to explain the correlations data nor is the concept of mini-jets necessarily 
required. 

In light of the NeXSPheRIO calculations, the highly structured correlation land- 
scape at RHIC should not necessarily be taken as an invalidation of the hydro- 
dynamic models. The correlations may simply reflect the need to abandon certain 
approximations, including the approximation of infinitely smooth initial conditions. 
Besides comparing to two-particle correlation data, these models can be used to cal- 
culate V2{2} and f2{4} to directly compare to data. It will be interesting to see how 
the correlation landscape in this model depends on the parameters of the model, in 
particular, the thermalization time and the freeze-out time. The connection of Vn 
fluctuations (related to two-parti cle c orrelations) to the lifetime of the system was 
first pointed out by Mishra et al \^^^\ In that reference the authors also introduce 
the anaology between yJv^ fluctuations and the power spectrum of the Cosmic 
Microwave Background Radiation. 

3.0.7. Addressing Uncertainties 

Uncertainties in the freeze-out prescription and the effects of the hadronic phase 
can be experimentally addressed through precise measurements of i^-mesons and 
f2-baryons. Models indicate that due to their small hadronic cross-sections, these 
hadrons are minimally influenced by the hadronic phase and reflect well the QGP 
phase. In addition, heavy flavor hadrons may help determine or provide a cross- 
check for the transport properties of the QGP. Another approach to extracting the 
viscosity is by studying the shape of V2{pt) versus system size. This approach does 
not rely on a model for the initial eccentricity. Uncertainties in the eccentricity 
and the initial conditions can be reduced through measurements of V2 fluctuations 
and two-particle correlations. These studies are ongoing. One can also measure u„ 
fluctuations for arbitrary n value. These are of course related to the two-particle 
correlation landscape which has already been extensively studied at RHIC. It will 
be of great interest to see how the correlation landscape predicted in hydrody- 
namic models with fluctuating initial conditions changes depending on the model 
parameters. The correlations data may help constrain quanties like the lifetime of 
the system. The studies listed above, along with a beam-energy scan at RHIC and 
the first data from LHC, will allow for more progress in understanding the matter 
created in heavy-ion collisions and its subsequent evolution. 

4. Summary 

In this review, V2 measurements were presented as a method for studying space- 
momentum correlations in heavy- ion collisions. The measurements of V2 indicate 
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Fig. 37. Hydrodynamic calculations including fluctuations in the initial conditions. The left panel 
shows the cvcnt-by-cvcnt V2 distribution. The right panel shows the two-particle correlations that 
arise due to the correlations in the initial energy density. 



the eccentricity in the initial overlap region is transferred efficiently to momentum- 
space. At top RHIC energy, the conversion is near that expectated from zero mean- 
free-path hydrodynamic predictions. The comparisons of data to hydrodynamics, 
however, depends on model calculations of the initial eccentricity. Several models 
for the initial eccentricity have been discussed. The mass, and Pt dependence of V2 
at < 1 GeV/c is found to be consistent with emission from a boosted source. 
Above that, the particle type dependence ofv2 exhibits a dependence on the number 
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of constituent quarks in the hadron, with baryons obtaining V2 values larger than 
mesons. 

The relationship between two-particle correlations, and V2 fluctuations has 
also been discussed. Calculations showing that some of the structures in two-particle 
correlations can be ascribed to fluctuations in the initial conditions, have been re- 
viewed. Measurements of correlations and V2 fluctuations can therefore be used 
to constrain models for the initial conditions. These constraints, along with im- 
proved measurements of the shape of V2{pt) as a function of system-size, improved 
measurements of and Q V2, measurements of V2 for heavy- flavor hadrons, mea- 
surements at LHC energies, and a beam-energy scan at RHIC will further improve 
our understanding of the properties of the matter created in heavy-ion collisions. 
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